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The variational cluster approach (VCA) based on the self-energy functional theory is applied to 
the two-dimensional symmetric periodic Anderson model at half filling. We calculate a variety of 
physical quantities including the staggered moments and single-particle spectra at zero temperature 
to show that the symmetry breaking due to antiferromagnetic ordering occurs in the strong coupling 
region, whereas in the weak coupling region, the Kondo insulating state without symmetry breaking 
is realized. The critical interaction strength is estimated. We thus demonstrate that the phase 
transition due to competition between antiferromagnetism and Kondo screening in the model can 
be described quantitatively by VCA. 
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I. INTRODUCTION 

The two-dimensional (2D) heavy fermion system has 
attracted much attention as a subsequent study of high- 
temperature superconductors and has recently been one 
of the central issues in the study of strongly correlated 
electron systems. For example, the heavy fermion mate- 
rial YbRh2Si2 shows a rapid change in the Hall coefficient 
as a function of magnetic field at zero temperature, which 
is accompanied by the antiferromagnetic (AF) to para- 
magnetic (PM) phase transition^ Also, heavy-fermion- 
like behavior is observed in the system of 3 He bi-layers 
adsorbed on graphite. 2 Generally speaking, competition 
between the magnetic ordering of localized spins through 
the Ruderman-Kittel-Kasuya-Yoshida (RKKY)? interac- 
tion and the nonmagnetic states induced by the Kondo 
screening 4 brings about the observed anomalous behav- 
iors in heavy fermion materials. From the theoretical 
point of view, the periodic Anderson model (PAM)£ is 
one of the simplified models for heavy fermion systems, 
which is believed to describe the competition between 
magnetic ordering and Kondo singlet formation observed 
in heavy fermion materials. 

In 2D quantum systems, the symmetry-broken mag- 
netically ordered state can be realized in the ground state 
at zero temperature and therefore one needs a method 
of calculation appropriate for infinite-size systems in the 
thermodynamic limit. In this paper, we therefore use the 
variational cluster approach (VCA)&£ based on the self- 
energy functional theory (SFT) 8 to consider the 2D pe- 
riodic Anderson model at half filling. Although the self- 
energies of the small clusters are used in the VCA calcu- 
lations and thus the long-range spin fluctuations beyond 
the cluster size are not taken into account, the quantum 
fluctuations within the cluster are treated exactly in this 
approach. We may therefore expect that this approach 
should be applicable to describe the possible symmetry 
breaking of the model in the thermodynamic limit be- 
yond the simple mean-field theory. We want to point 



out that the present calculation is the first one where the 
VCA is applied to the 2D PAM, as far as we know. 

We will show that, by means of VCA, the symmetry 
breaking due to the AF ordering of localized spins oc- 
curs in the strong coupling region, whereas in the weak 
coupling region, the Kondo insulator without symmetry 
breaking is realized. The critical interaction strength 
will thereby be determined. We will also calculate the 
staggered magnetic moment as a function of the inter- 
action strength and show that the phase transition is 
of the second order. We will furthermore calculate the 
single-particle spectra and densities of states (DOS) to 
discuss the effects of electron correlation on the quasipar- 
ticle band structure. We will thus show how the system 
changes from the AF insulator, Kondo insulator, to the 
band insulator, with decreasing the interaction strength. 

This paper is organized as follows. In Sec. [Ill we 
present our model and method of calculation. In Sec lIIIi 
we present our results of calculations for the stability 
of the AF ordering, staggered magnetic moment, single- 
particle spectra, and DOS by VCA. We summarize our 
work in Sec. [IVJ 



II. MODEL AND METHOD 
A. Model 

We consider the PAM defined on the 2D square lattice. 
The Hamiltonian is given by 

H = -t 5>.W + H.c.) - Vj2(4ji° + H.c.) 

(ij) i<7 
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where c 4(T (fia) is the annihilation operator of an elec- 
tron at site i with spin a in the conduction-electron c 
(/-electron /) orbital, and n{ a = f} a ficr is the electron 
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number operator in the / orbital, t is the hopping pa- 
rameter between the nearest-neighbor c orbitals, V is the 
on-site hybridization parameter between the c and / or- 
bitals, U is the on-site repulsion on the / orbital, and ej 
is the energy level of the / orbital with respect to that 
of the c orbital set to be the origin of energy. In the fol- 
lowing calculations, we consider the symmetric case, i.e., 
the case with £/ = —U/2. We also focus on the electron 
densities at half filling, i.e., 2N S electrons in the N a unit 
cells, where the unit cell contains one c and one / orbital. 
We hereafter set t = V = 1 as the unit of energy and we 
change the value of the interaction strength U. 

B. Variational cluster approach 

Let us first briefly review the formulation of SFT 8 and 
present the method of calculation of the magnetic order- 
ing by VCA&i in order to make our paper self-contained. 

We consider the system of the Hamiltonian H = 
Ho(t) + Hi(U), where t and U denote the one-particle 
and interaction parameters of H, respectively. In general, 
the grand potential is given from the stationary point of 
the self-energy functional 

n[S]=F[S]+Trln[-(Go 1 -S)- 1 ], (2) 

where F[T] and Go are the Legendre transform of the 
Luttinger-Ward potential $[G] and the bare Green func- 
tion, respectively. The rigorous variational principle 
5f2[£]/<5£ = gives the Dyson equation G -1 = G^ 1 - S, 
where G is the physical Green function. 

In the above expression (2), F[T] is a universal func- 
tional of the self-energy; i.e., F[S\ remains unchanged 
for an arbitrary reference system of the Hamiltonian 
H' = Ho{t') + Hi(U) that has the same interaction 
part as the original system has, but with modified one- 
particle parameters. We here introduce the restriction of 
the space of the exact self-energies of the original system 
to the set of exact self-energies of the reference system. 
Because of this restriction, the following procedure be- 
comes approximate but it enables us to obtain the grand 
potential of the original system from the stationary point 
of the £(£') functional 

r>[s(t')] = n' + TrH-iGo 1 - e(0) _1 ] 

- TrlnHG'o 1 -^'))- 1 ], (3) 

where Q', S(t'), and G'o are the grand potential, exact 
self-energy, and bare Green function of the reference sys- 
tem, respectively. The condition dfl[Y,(t')]/dt' = gives 
an appropriate reference system that describes the origi- 
nal system approximately. 

In VCA, we first divide the original infinite lattice into 
the finite-size identical clusters. By switching off the hop- 
ping parameters between the identical clusters, we con- 
struct the reference system as an assembly of the exactly 
solvable finite-site clusters. One of the major advantages 
of VCA is its ability to describe the symmetry-breaking 



long-range order by introducing suitably chosen fictitious 
symmetry-breaking Weiss fields in the set of variational 
parameters t' . In order to discuss the competition be- 
tween the AF ordering and Kondo screening in the pa- 
rameter space, we here introduce staggered magnetic field 
h' on the / orbitals in the cluster Hamiltonian as a vari- 
ational parameter. We thus obtain the Hamiltonian of 
the reference system, H 7 , which is given by 

H' = J2 H '*> W 

R 

h' r = -*5>U* + H - c -) - v T,(4ji« + h.c.) 

(ij) icr 
i icr 

+h 'Y^e i ^(n{ 1 -n f il ), (5) 

i 

where R is the label of the clusters, i and j are the labels 
of the sites within the cluster R, and Q = (ir,ir). 

In the present study, we use a 6-site (2 x 3) cluster to 
search for the stationary point of 0[E(/i')] with a condi- 
tion dQ[H(h')]/dh' — as discussed above. We should 
note that the shape of the cluster introduced as a refer- 
ence system is not commensurate with the AF ordering. 
We therefore treat a 12-site (2 x 6) cluster as a supercell 
by combining the two 6-site clusters. We treat the inter- 
cluster hopping elements as well as the hopping elements 
between the supercells "perturbatively" p i.e., we use the 
self-energies of the 6-site clusters to calculate the Green 
function of the original infinite system as well as that of 
the reference systems (an assembly of the identical 12-site 
clusters) via the Dyson equation and obtain the values 
of for various values of h! by using the Eq. ([3]). 

III. RESULTS OF CALCULATIONS 
A. Stability of the antiferromagnetic ordering 

It is known that the AF ordering of the / electrons is 
realized in the ground state of the strong coupling region 
of PAM (as well as Kondo-lattice model) in 2B^lLi2. We 
demonstrate this in Fig. [TJ where the calculated values of 
f2[£(/i')] — fi[£(0)] per site for several values of U near the 
critical point are shown. We find the following: (i) The 
value has a minimum at a finite value of h' for U > t/ cr , 
which indicates that the symmetry-broken AF ordering 
is stabilized for U > U cx . (ii) The value of h! at which 
f2[E(/i')] — f2[S(0)] has a minimum approaches with 
decreasing U to U — » U cr . (iii) The critical value of U 
is determined as U cr — 2.7. This value is comparable to 
(but is slightly smaller than) the result of the quantum 
Monte Carlo calculation 10 where the value U CI ~ 2.95 
is reported. The reason of the overestimation of the AF 
stability in VCA may be explained as follows: In the VCA 
calculation, we use the cluster representable self-energies, 
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FIG. 1: Calculated results for fi[£(/i')] - fi[£(0)] (per site). 
The right panel shows an enlargement of the small h' region 
of the left panel. We show the results for several values of U 
near the phase transition. Dotted horizontal line is a guide 
for eyes. 



i.e., exact self-energies of small clusters, as the trial self- 
energies. Thus, the long wave-length spin fluctuations 
beyond the cluster size are not taken into account^ Also, 
we use the staggered magnetic field on the / orbitals h' 
as a single variational parameter, which suppresses the 
quantum spin fluctuations. This may also be responsible 
for the overestimation. (iv) For U < U CT , we find the 
minimum at hi = 0, which indicates that there is no 
long-range AF ordering in the system. 
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where f2 is the grand potential of the system per site and 
^cxt (^ext) i s the external staggered magnetic field act- 
ing on the / (c) orbitals. We find the following: (i) The 
staggered magnetic moment of the / electrons decreases 
continuously to when we decrease the value of U from 
the strong coupling region to U CI . Thus, the phase tran- 
sition is of the second order, (ii) The staggered moment 
of the c electrons also shows the similar behavior but the 
polarization is of the opposite sign, (iii) The obtained 
staggered moments may be overestimated as in the case 
of the Hubbard model in 2D, 7 which is due again to the 
overestimation of the stability of the AF ordered state in 
VCA as discussed above. 

We also point out that the calculated staggered mo- 
ments (Mf) and (m c ) are found to show a power-low 
behavior in the vicinity of the transition point as ~ 
(U—U CV )P with the exponent (3 ~ 0.5, which is consistent 
with the mean-field value (3 — 0.5 within the numerical 
accuracy. Thus, the VCA calculation, which neglects the 
long-range spin fluctuations beyond the cluster size, gives 
the results equivalent to those of the simple mean-field 
approximation at least in the description of the critical 
behaviors such as the critical exponent. 



B. Staggered magnetic moment 
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FIG. 2: Calculated results for the staggered magnetic moment 
of the / orbitals (Mf) and c orbitals (m c ) as a function of U. 
Dotted line represents the critical value f/ C r = 2.7. 



In Fig. we show calculated results for the staggered 
magnetic moment of the / orbitals (Mf) and c orbitals 



C. Single-particle spectra and densities of states 

We then calculate the single-particle spectra 13 defined 
by the imaginary part of the "Fourier transform" of the 
optimized physical Green function^ We also calculate the 
DOS from the fe-space integration of the imaginary part 
of the optimized physical Green function. The results 
are shown in Fig. [3J where we use the AF Brillouin zone 
when the system is in the symmetry-broken AF state (see 
Figs. [3j(a) and (b)), but for nonmagnetic states, we use 
the standard first Brillouin zone (see Fig. [3J(c)). 

In Fig. [3J (a), i.e., for U = 10, we can first identify 
the "upper and lower Hubbard bands" for the / elec- 
trons, which are almost dispersionless and are separated 
by an energy ~{7. We can also identify the lower-energy 
dispersive bands in Fig. [3] (a). Here, we use the spin- 
density-wave (SDW) dispersion to fit the spectra. The 
SDW dispersion can be obtained by diagonalizing the 
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FIG. 3: (Color online) Calculated results for the single-particle spectra (left and middle panels) and DOS (right panels) for 
(a) U — 10, (b) (7 = 4, and (c) U — 2, where uj = corresponds to the Fermi energy. The left and middle panels show the 
spectra of the c and / electrons, respectively. In the right panel, the solid and dotted curves are the DOS for the / and c 
orbitals, respectively. The artificial Lorentzian broadening of r\ = 0.05 is included. 



SDW Hamiltonian T/sdw defined by 



Hsdw = 
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where A and B are the sublattice indices, V and Ef are 
the effective hybridization parameter and effective en- 
ergy level of the / orbital, respectively, Mf (m c ) is the 
staggered magnetic moment of the / (c) orbitals, and 
£fc = —2t(cosk x + cos k y ). We assume Mf and m c to 
have the values obtained in Eqs. (6a) and (6b) and we 
fix Ef to be 0. We determine the value of V so as to re- 
produce the size of the SDW gap. We find that the fitting 
works well for the dispersions of the lower-energy bands 
but the spectral weight on the / orbital differs very much 
from that of the VCA calculations since the upper and 
lower Hubbard bands for the / electrons do not appear 
in the SDW spectral functions. We then find the value 



V ~ 0.35 from the fitting, indicating that the quasipar- 
ticle is not quite heavy. In other words, with increasing 
U, the AF ordering occurs in 2D before the quasiparticle 
mass is strongly enhanced. 

In Fig. [3](b), i.e., for U = 4, we find that the localized 
energy level is not well defined but there is a band repul- 
sion in the spectra at ej = ±17/2. The spectral weight 
of the / electrons becomes large near the Fermi energy 
for all the momenta. Also, by comparing the results of 
the non- interacting case (U — 0), the sharp peak just 
below the Fermi energy is observed in the partial DOS 
of the / orbital (see the right panel of Fig. [3] (b)). Thus, 
we conclude that this peak not only arises from the hy- 
bridization but is caused by the many-body resonance, 
which corresponds to the Kondo resonance peak in the 
metallic state. 

In Fig. [3] (c), i.e., for U = 2, where there is no AF 
ordering in the system, we find that the spectra look very 
similar to the spectra of non-interacting case. However, 
we again find that the localized energy level is not well 
defined but there is a weak band repulsion in the spectra 
at ef = ±17/2. 
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D. Charge gap and Spin gap 
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FIG. 4: Calculated result for the ratio of the spin gap to 
charge gap A s /A c as a function of U (U < U C r)- The 8-site 
cluster is used. 



To clarify the behavior in the weak-coupling region 
where there is no AF ordering, i.e., U < U a , we calcu- 
late the spin and charge gaps defined as A s = Eq(N^ + 
1,N L - 1) - Eq (JVf , N± ) and A c - [£ (A T + l.JVj. + 
1) + Eq(Ni - l,iV; - 1) - 2E (N h N l )]/2, respectively, 
where E (N^,Ni) is the ground-state energy of a cluster 
with iVf up-spin and JVj_ down-spin electrons. Because 
the two-particle Green functions cannot be calculated di- 
rectly from VCA, we here use an exact-diagonalization 
technique on small clusters. We use the 8-site, 16-orbital 
cluster with periodic boundary condition to calculate the 
ground-state energies and estimate the spin and charge 
gaps. In Fig. QJ we show the ratio of the spin gap to 
the charge gap A s /A c thus obtained as a function of U, 
where the result only at U < U cr (with U CI determined in 
Sec. Ill A) is shown because no phase transition occurs 
in finite-size systems. We find A c > A s for all values 
of U (< Ucr)) indicating the system to be in the regime 
of the Kondo insulator ji£ i.e., there is no long-range AF 
ordering, where localized spins are screened by the for- 
mation of the Kondo singlet state. As U decreases to 
0, we find that the two gaps tend smoothly to the same 
value, i.e., A s /A c — » 1, indicating the system tends to 
the non-interacting band insulator. 
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IV. SUMMARY 

In summary, we have applied the VCA based on the 
SFT for the first time to consider the symmetric PAM at 
half- filling in 2D. We have thus demonstrated the validity 
of the approach by discussing in particular the competi- 
tion between antiferromagnetism and Kondo screening 
in the thermodynamic limit at zero temperature. We 
have shown that the symmetry-broken AF ordering of 
localized spins is realized in the strong coupling region 
U > U CI and the Kondo insulating behavior is realized in 
the weak coupling region U < U CI . We have determined 
the critical interaction strength as U CI = 2.7. We have 
calculated the staggered magnetic moment as a function 
of the interaction strength and have shown that the phase 
transition is of the second order. We have also calculated 
the single-particle spectra and density of states. We have 
thereby discussed the effect of electron correlations on the 
quasiparticle band structure. We have applied an exact- 
diagonalization technique on small clusters to calculate 
the ratio of the spin gap to charge gap in the weak cou- 
pling region and found that the Kondo insulating state 
continuously tends to the non-interacting band insulator 
with decreasing the value of U to 0. 

We thus have shown that the present approach is very 
useful for considering the electronic states of PAM in 2D. 
To improve the accuracy of our results, one may intro- 
duce additional variational parameters, such as the hop- 
ping terms, to take into account the quantum fluctua- 
tions more efficiently and suppress the overestimation of 
the stability of the AF ordering, which we want to leave 
for future studies. 
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